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4  &  Space  Propulsion  and  Power  Portfolio 

Research  to  Understand,  Predict,  and  Control  Complex 
interactions  of  the  Matters  in  Space  Propulsion  Systems 
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Micro  plasmas 


Coupled  Materials  and 
Plasma  Processes  Far  From 
Equilibrium 


magnetron  Sandia  Saturn  Pulsed 
Power  Generator 


Satellite  contamination 

■  HI 


DUAL-  MODE  PROPULSION  -  micro  chemical  thruster 


Electrospray  propulsion 


Novel  Energetic  Materials 

Reduced  Basis  and  Stochastic  Modeling 
of  high  pressure  combustion  dynamics 
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Space  Propulsion  and  Power  Portfolio 

Program  Interactions  and  Trends 
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Coupled  Materials  and  Plasma  Processes 
Far  From  Equilibrium-Material  /  Plasma 

Coupling 


Sheath  Potential 
Instabilities 
~4  MHz 


High  Frequency  Bulk 
Plasma  Oscillations 


Coherent  structures 


Inverse 
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width  of  the  power  spectra  in  frequency  space 
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Wall  Material 
(SEE  emitter) 
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50  60x10 

Bulk  Plasma 


Plasma  Sheath 

Why  do  we  care? 


Coherent 
instability 
in  Hall 
Thruster 


o  Change  the  particle  &  energy  fluxes  ! 
o  Change  the  stability  and  lifetime  of  device  ! 
o  Change  the  characteristics  (SSA)  of  device  ! 
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Y  Y  Wall  Material  Experiments  (Raitses,  Princeton) 

Wall  Architectures  Have  Significant  Effects  on  Secondary 
Electron  Emission  (SEE)  and  Discharge  Behavior 


Scanning  Electron  Microscope  Images 

“Black”  regions  -  No  SEE:  “White”  regions  -  SEE 


Velvet  fibers  on  the  surface 
strong  SEE  from  fibers 


Velvet  fibers  perpendicular  to 
the  surface-No  SEE  except 
from  fiber  tips 
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•Need  to  take  into  account  spatial  /temporal  variations  of  plasma  scale,  XD  (Debye  length),  due  to  plasma 
design  and  instabilities 


•  Y.  Raitses,  I.  D.  Kaganovich,  A.  Khrabrov,  D.  Sydorenko,  N.  J.  Fisch,  A.  Smolyakov,  IEEE  Transactions  on  Plasma  Science  39,  995  (2011  ) 
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Modeling  of  plasma-wall  interaction  (Kaganovich.  Princeton) 

^St^Kinetic  simulations  revealed  a  new  regime  of  plasma-wall 
interaction  with  a  very  strong  secondary  electron  emission 


Beams  of  SEE  electrons 


Left  wall 


SEE  yield  is  the  number  of 
secondary  electrons  emitted 
per  incident  primary  electron 


Y  1 


L _ V. 


Energy  of  incident  electron,  eV 

Effective  secondary 

* 

electron  emission  y 

accounts  for  non- 
Maxwell  ian  effects 


Experimentally  validated  by 
N.CIaire  and  F.  Doveil,  Aix- 
Marseille  Universite/CNRS 


t(10'6s) 


•  electrons  acquire  enough  energy  from  the  electric  field  parallel 
to  the  wall,  causing  ,  y*  >  1 

•  Sheath  collapse  leads  to  extreme  wall  heating  by  plasma  and 
plasma  losses  (bad) 


M.  D.  Campanell,  A.V.  Khrabrov,  I.  D.  Kaganovich,  Physics  of  Plasmas  19, 123513  (2012) 
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^  k  Small-Scale  Turbulence  and  Inverse  Cascade 

Generating  Large  Scale  Coherent  Structures  for  Hall 
discharges  (magnetron  similar) 


■  Generated  by  nonlinear  coupling  between  fast-growing  unstable  higher  frequency  modes  and  lower  frequency  modes 

■  These  large-scale  fluctuations  correlate  to  measured  transport  properties 


width  of  the  power  spectra  in  frequency  space 
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A  A  Cappelli  /  (Stanford).  Haraus  (AFRL/RQRS) 

Successfully  synchronized  continuous  wave  laser 
induced  fluorescence  to  coherent  structures  in  plasmas 

•  N.  A.  MacDonald,  M.  A.  Cappelli,  and  W.  A.  Hargus,  Rev.  Sci.  Instruments  83, 113506  (2012) 

•  Time-Synchronized  continuous 
wave  LIF  Measurements  of  Velocity 

-  First  Results 
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Velocity  fluctuations  vastly  different 
depending  on  probed  location  in  magnetized 
plasma 

Physical  probes  can  only  measure  the 
average  velocity! 
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Wirz  /  UCLA-  YIP-  Center  of  Excellence  (Univ.  of  Michigan) 

Self-consistent  analytical  theory  used  to 
optimize  micro  plasma  source 


*  Critical  research  for  in-space  propulsion,  plasma-enhanced  combustion,  plasma  aerodynamics,  small  scale  sensors, 
directed  energy  devices,  plasma  processing  ... 


Goal:  Develop  efficient  and 
stable  cusp-confined  micro 
discharge 

•  New  “stream  function” 
magnetic  field  analysis 
developed  to  predict  plasma 
behavior  throughout  discharge 


Optimization  using  analytical  theory  yields  2X  increase  in  plasma  production  and  stable 

plasma  over  baseline 

•  Model  trends  experimentally  validated 
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•  Mao,  H-S.,  Wirz  R.  E.,  Appl.  Phys.  Lett,  accepted  Nov  2012 
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4  k  Electrospravs.  dual-mode  propulsion  (Lozano-MIT-YIP) 

^^Molecular  Dynamics  Simulations  reveal  how 
to  mitigate  Fragmentation  of  Solvated  Ions 


When  electric  field  (or  electric  pressure)  >  surface  tension  +  internal  liquid  pressure 
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Single  Porous  emitter 


Full-atom  MD  with  E  field 


•  MD  results  suggest  the  use  of 
complex  ions,  as  there  are  more 
degrees  of  freedom  to  dissipate 
internal  energy 


EMI-BF4 


Less  complex  dimer 
100%  fragmentation 


pmi-pf6 


More  complex  dimer 
55%  fragmentation 


Carla  S.  Perez-Martineza)  and  Paulo  C.  Lozano,  J.  Vac.  Sci.  Technol.  B 
30(6),  Nov/Dec  2012 

■  Ofetrityution  Ar  Ttpprwect  forpublic  release;  distribution  is  unlimited 


4  4  Center  of  Excellence  (Univ.  of  Michigan)  -  King  (Michigan  Tech) 

Can  we  obtain  Taylor  cone  without  a  physical 

emitter  or  capillary? 


Taylor  cone  forming  at  tip 
of  physical  emitter 


Electric  Field 
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Tension 


Physical  needle  enhances  E-field  - 
Taylor  cone  forms. 


Taylor  cone  forming  at  the  peak 
of  Normal  Field  Instability 


Surface 


Fe203 

nanoparticle 


Permanent  magnet 


Concentrated  B  at  crest  attracts  more 
ferromagnetic  fluid  to  crest,  amplifying  the 
perturbation  into  an  instability.  Crest  enhances 
E-field  -  Taylor  cone  forms. 


Ionic  liquid  becomes  superparamagnetic  with  an  addition  of 


ferromagnetic  Fe203  nanoparticles! 
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Center  of  Excellence  (Univ.  of  Michigan)  -  King  (Michigan  Tech) 

First  stable  ionic-liquid  ferrofluid  synthesized  by 
Hawkette,  et  al  2010  (Australia) 


•  Jain,  N,  Wang,  Y,  Jones,  SK,  Hawkett,  BS,  and  Warr,  GG,  “ Optimized  steric  stabilization  of  aqueous  ferrofluids  and  magnetic  nanoparticles,  "Langmuir,  26(6),  2010,  pp.  4465-4472. 
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^  1  Center  of  Excellence  (Univ.  of  Michigan)  -  King  (Michigan  Tech) 

^4  Progress  -  Normal  field  Instability  with  E-field 


"Normal  field  instability"  peaks  (no  electric  field)  can  be  When  an  electric  field  is  applied  these  peaks  will  get 

rounded  or  they  can  be  very  sharp  (shown  in  inset)  even  sharper 


•  Weak  magnetic  field  with  high-surface-tension  ionic  liquid  ferrofluid  shown.  Normal- 
field  instability  peaks  (left)  are  rounded  in  this  case. 

•  Application  of  moderate  electric  field  sharpens  peaks  due  to  electrostatic  force  (right) 
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Center  of  Excellence  (Univ.  of  Michigan)  -  King  (Michigan  Tech) 

Challenges 


Q1 :  Emitter  tip  density  (tips  per  mm2) 
is  a  function  of  surface  tension, 
nanoparticle  magnetization,  and 
magnetic  field.  Can  we  synthesize 
an  “ideal”  ILFF  for  space  propulsion 
that  has  high  thrust  density  (milli- 
Newtons  per  mm2)? 


Q2:  An  ILFF  is  a  complex 
ferromagnetic,  electrically  conductive, 
and  non-homogeneous  colloidal  fluid. 
What  are  the  velocities  of  the 
molecular  and  macro  species  emitted 
when  the  fluid  is  electrosprayed  and 
what  is  the  beam  divergence? 


Mass  Spec  of  201 20830  Emim-Ac  (SS)  Cation  2800V  20deg 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited 


w 


Space  Propulsion  and  Power 


A  Brief  History  of  Nanoenergetic  Material 


Examples  from  AFOSR  program: 

•  1st  Generation 

■  Nanometer-sized  Al  powder/conventional 
propellants 

■  Some  performance  gain,  variable  results 


•  2nd  Generation  -  7op  down  approach 

■  Quasi-ordered  nanometer-sized  inclusions  in 
energetic  matrix 

■  Coated  nanometer-sized  metal  powders 

■  Controlled  oxidation,  improved  storage  lifetime 


AI  + 

fluoropolymer  Propellant  w/ Baseline  Al 


Modified 

Al 


Nano-AI  Encapsulated  with 
Ammonium  Perchlorate 


•  3rd  Generation  -  bottom-up  approach 


Organized  multiscale  processing  to  enable  the  insertion  of  nanoenergetic 
materials  into  larger  units 

3-dimensional  nanoenergetics  for  reaction  control 
Controlled  reactivity,  Improved  manufacturability/processing 
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2nd  Generation: 

Aluminum  with  nanoscale  inclusions  of  fluoropolymers 
improves  propellant  performance  (Son/Purdue) 


2012 


Al  +  fluoropolymer 
Powder  Mixture 


Uniform  composition  and 
induced  lattice  defects  (stored 
energy) 


•  Experiments  validate  the  hypothesis  that  aluminum  with  nanoscale 
fluorocarbon  inclusions  in  a  solid  propellant  yield  reduced  ignition 
temperature,  increased  burning  rates,  and  smaller  agglomerate  size 

2  f  : 

□  Baseline,  (15-75 /inn) 


Propellant  w/  Baseline  Al  w /  Modified  Al 


•  Other  nano-inclusion  materials  will  be  explored,  including  piezoelectric  polymers  to  achieve 
smart/functional  control  of  rate  or  sensitivity 


•  Travis  R.  Sippel,  Steven  F.  Son,  and  Lori  J.  Groven,  “Altering  Reactivity  of  Aluminum  with  Selective  Inclusion  of  Polytetrafluoroethylene 
through  Mechanical  Activation”  Propellants,  Explosives,  Pyrotechnics,  2013 
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1  t  3rd  Generation: 

^^^yFY  2012  MURI  :  Smart,  Functional  Nanoenergetics  Design  from 
the  atomistic  /  molecular  scale  through  the  mesoscale 


biological  example 

of  multiscale  effects 


Gecko 


Rows  of  setae  from 
a  toe  (micron) 


Spatulae  (nm) 


i=>  Multiscale  Energetic  Composites  Fabricated  on  Porous 

Silicone  Substrates  (Yetter/  Penn  State) 


Autumn,  K.,  eta!.,  Nature, 
405,  681-684,  2000 


Side  View  -  Reaction  Propagation 


pillars  were  ~  35  pm 
tall  and  have  8  pm 
square  bases 
separated  by  ~8  pm 


Pore  Diameters  ~  20nm 
and  filled  with  oxidizer 

Mg(CI04)2 


•  V.K.  Parimi,  S.A.  Tadigadapa,  and  R.A.  Yetter,  J.  Micromechanics  and  Microengineering,  22,  5,  2012 
Distribution  A:  Approved  for  public  release;  distribution  is  unlimited 


Third  Generation: 

Organized  Multiscale  Energetic  Al  Composites  for 

Combustion  Control 


Step  1 : 

Generation 

of 

passivated 
Al  clusters 


AI-77  Cluster 


E 
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Pyrophoric  and 
extremely  fast 
reacting  Al  particle 

- > 


Metalloid  Al  clusters  modeled  after  Schnockel 
aluminum  cluster  chemistry  (Ecker,  A.,  Weckert,  E., 
Schnockel,  H.,  Nature  1997,  387,  379) 


Nanoscale  passivated  Al  Cluster  using  ligands 


Bond  energy  at  ligand 
head  controls  IGNITION 


Passivating  ligand 

provide  additional  energy 
during  combustion 


Al  cluster  TEMPERATURE  of  cluster 


Nanoscale 
passivated  Al 
clusters  to  be 
evaluated  as 
additives  to 
liquid  fuels  <& 
propellants 


Aerosol  of  Al  cluster  and  gas 

generator  molecules 


Step  2: 

Generation 

of 

mesoscale 

energetic 


composites 


Energetic  gas 
generators  with  large 
AH  of  combustion 


Mesoscale  composite  of  Al  cluster  and  gas  generator 

for  incorporation  in  solid  and  liquid  propellants 


Aerosol  based  self 
assembly-  a  bulk 
manufacturing  process 
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Gasification  and  decomposition  temperatures 
— >  gas  generator  molecules  control  s  energy  and 
temperature  of  mesoparticle  disruption 

Distribution  A:  Approved  for  public  release;  distribution  is  unlimited 


“saa0 

J 


disassemble  and 
release  highly 
dispersed  reactive 
nanostructures  at  a 
predetermined 
pressure  and 
temperature  for 

controlled 

combustion 
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Summary  and  New  Research  Areas 


Coupled  Materials  and  Plasma 
Processes  Far  From  Equilibrium 


Control  of  Coherent  Structures: 

•  Stable  Plasma  Propulsion  with 
extended  lifetime 

•  Plasma  Photonic  Crystals 


Electrosprays 

x 

Novel  Energetic  Materials 


•  Non-Equilibrium  Plasmas  in  Liquid 
Propellants 

•  Miscible  liquid  boranes  in  propellants 

•  Transformable  Energetic  Materials  via 
co-crystallization  or  molecular  blending 


•  Quantum  Lattice 
Modeling  for 

Reduced  Basis  and  Stochastic  Modeling  of  high  Multiphase  reacting 

pressure  combustion  dynamics  flows 
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Space  Propulsion  and  Power  -  New  Research  Areas 
Control  of  the  Coherent  Periodic  Microstructures  in  plasmas 
may  lead  to  reconfigurable  THz  Plasma  Photonic  Crystals 


From  micro  plasmas 


Dielectric 


Electrode 


Electromagnetic 

wave 

- > 
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Dielectric  Electrode 


Discharge  pulse 


Advantages  of  using  Plasma  Photonic  Crystals  : 

■  Tunabilitv  :  variable  Refractive  Index  and  Band  Gap 

■  Reconfigurable  structure  :  variable  crystal  Geometry/Symmetry 

■  Inherent  Nonlinearity  :  Harmonic  Wave  Generation 


From  bulk  coherent  structures 


CHALLENGES: 


Can  researchers  achieve  plasma  crystal  lattice  scale  to  THz 
electromagnetic  wavelength  with  a  plasma  frequency  close  to 
electromagnetic  frequency  ? 

Can  researchers  generate  and  stabilize  coherent  organized 
microstructures  from  large  otherwise  uniform  dense  plasmas  ? 
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Space  Propulsion  and  Power  -  New  Research  Areas 

Non-equilibrium  plasmas  in  liquid  propellants 


Recently,  non-equilibrium  plasmas  have  been  formed  in  liquid  water  without  formation  of 
gas  bubbles  yielding  propagation  velocities  of  5000km/s  for  low  energy  nanosecond 
discharges 


(Starikovskiy  et  al.,  Plasma  Sources  Sci.  and  Technol.  20, 1,  2011). 


wm 

o 

o 

t=  250  ps 

t  =  500  ps 

t  =  7  SO  p$ 

t  =  1000  ps 

Impact: 

The  formation  of  highly  ionized  channels  in  condensed  media  without  void  formation 
may  create  “liquid  plasma”  applications  for  ignition  assistance,  in-situ  propellant 
modification,  and  accelerated  combustion. 


Challenge: 

However ;  the  roles  of  low  energy  non-equilibrium  discharges  with  propellants 
and  models  for  describing  discharges  in  dense  media  are  poorly  understood 
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w  Space  Propulsion  and  Power  -  New  Research  Areas 

Miscible  liquid  Boranes  in  RP1  (kerosene) 


Borane-N,N  diisopropyl 
ethylamine  complex 


Borane-triethylamine  complex 


Impact:  low-toxicity  and  combustion  instability  control 

Challenges: 

•  A  fundamental  study  is  needed  to  synthesize  optimal  amine-boranes  and/or 
other  borane  molecules  that  are  miscible  in  RP-1  with  coupled  simulation 
and  feedback  from  small-scale  characterization  to  develop  optimal 
propellants  and  additives 

•  Applicable  to  ionic  liquids? 
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Space  Propulsion  and  Power  -  New  Research  Areas 

Transformable  Energetic  Materials 


Example:  Co-Crystallization 


Material  A 


Material  B 


*  *  * 


high  sensitivity  low  density 

good  poor  performance 

performance 


Insensitive  cocrvstal  of  CL-20  and  TNT 


co-crystallization 


low  sensitivity 
good  performance 
propellant 

ability  to  transform  when 
stimulated 


- > 

Application 

of 

“switching” 
energy 
(electro¬ 
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Bolton,  O.;  Matzger,  A.J.  Angew.  Chemie  Int.  Ed.  2011,  50,  8960-8963 


Challenges: 


•  Can  we  engineer  green  energetic  materials  that  can  transform  (e.g. 
from  propellant  to  explosive,  or  even  in  situ  modification)? 


•  Can  we  model  and  predict  transformations  and  necessary  stimuli? 
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Space  Propulsion  and  Power  -  New  Research  Areas 

Quantum  Lattice  Algorithms  for  the  Simulation  of 
Multiphase/Multicomponent  Fluid  Phenomena 


Classical  computer  bit: 
Quantum  qubit: 
Quantum  advantage: 


0  or  1 

|#)=  |0)  +  1 1)  (quantum  state) 

instantaneous  result,  noiseless,  unconditionally 
stable,  reversible 


One  or  more  qubits  can  be  arranged  on  lattice  sites  (here  2-D) 


^ 
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o.  ,tx  : 


System  state  at  time  t 


System  evolves  in  time  by 

|  (xv...,xn\t  +  ))  =  SC|  (xv...,xnU)) 

where  S  is  a  streaming  operator  (flow) 

q  is  a  collision  operator  (fluid  interactions) 
both  implemented  via  quantum  gates 


Algorithm  can  be  run  on  quantum  or  classical  computers 

Shown  by  Yepez1  to  reproduce  the  lattice  Boltzmann 
equation,  a  popular  CFD  methodology. 


^epez,  J.,  “Quantum  Lattice-Gas  Model  for  Computational  Fluid 
Dynamics,”  Physical  Review  E,  Vol.  63,  2001,  p.  046702. 
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Plasma  Metamaterials  and 
Plasma  Optics 


plasma  metamaterials,  we  can  use  the  effects  of  Imfs) 

Refractive  index  N  is  imaginary  - 
Bulk  Electromagnetic  waves  cannot  propogate 
But  surface  plasmons  possible 


Refractive  index  N  >  0 
Low  density  collisionless 
Plasmas  -  photonic  crystals  with 
tunable  band  gaps  possible 


Real  part  of 
emissivity 


http://en.wikipedia.org/wiki/Ne 

gative_index_metamaterials 


N  is  negative  - 
Metamaterial  causes  light  to 
refract,  or  bend,  differently 
than  in  more  common 
positive  refractive  index 
materials 


N  is  imaginary- 
Gyrotropic  material 
(with  imposed  magnetic 
fields),  leading  to 
Faraday  rotation  and 
Optical  Kerr  effect,  one- 


light 

source 


polarizer 


way  waveguides 
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Plasma  Photonic  Crystals 


A  plasma  photonic  crystal  (PC)  is  an 
array  of  plasma  structures  which  have 
the  unique  ability  to  control  the 
propagation  of  EM  waves 


Novel  EM  response  can  be  obtained 
such  as  wave  guiding,  and  spectral 
filtering 


EM  Waves 

JW'- 

JWV 


Regulation  (on/off  state, 
parameters} 

Transmitted  waves 

*  *  * 

J  *  *  *  * 


Periodic  Micro  plasmas 


O.  Sakai  and  K.  Tachibana  (2012) 


Dielectric  Electrode 


Dielectric  Electrode 
O.  Sakai  and  K.  Tachibana  (2012) 
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Plasma  Photonic  Crystals 


State-of-the  art  PPC’s  are  designed 
for  controlling  several  mm-wavelength 
radiation  (tens  of  GHz) 

■  Plasma  scales  -  several  mm’s 

■  Plasma  densities  ~1 013  cnr3 

Propagation  bands  and  bandgaps 
appear  near  cut-offs  and  resonances 
(e.g.,  plasma  frequency) 

Future  trends  and  applications  seek 
to  control  or  manipulate  higher 
frequencies  -  hence  higher  plasma 
densities  and  smaller  plasma  scales 


1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0.0  0.1  0.2  0.3  0.4  0.5 


n  1  r 


plasma  vacuum 

\  \ 


Band  Gap 


Band  Gap 


Imaginary  Wavenumber 
Real  Wavenumber 


-i - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1- 


kal  2  7i 

Normal  incidence  dispersion  in  a  ID  plasma-air  photonic 
crystal.  The  inset  shows  the  configuration.  The  plasma 
density  is  ne  =  1014  cm~3,  dp  =  1.20  mm  and  lattice 
parameter ,  a  =  3  mm. 


Terahertz  EM  Wave  Manipulation 


ne  sl015-10"W3 
a  «  d  ~  0.1  mm 
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Plasma  Photonic  Crystals 


Externally  generating  an  array  of  microplasmas 
(using  hollow  anode  discharge  arrays,  for  example) 
will  become  increasingly  difficult,  particularly  for  the 
smaller  scales  and  higher  plasma  densities  needed 
for  controlling  THz  and  FIR  waves 

(a) 

Alumina 


Anode 


Cathode 


•  e  i 

e  €> 

•  e  < 

o  • 

©  •  G  €> 

g  e  © 

•  •  ©  0 


Martin  et  al  (2012) 

Plasmas  (and  magnetized  plasmas  in 
particular)  constitute  highly  non-linear  and 
naturally  unstable  systems. 


Seemingly  uniform  plasmas  (see  top  left  figure 
to  the  left)  can  be  “coaxed”  into  self-organizing 
into  regular  patterns 


Challenge: 

Can  researchers  generate  and  stabilize 
coherent  organized  microstructures  from  large 
otherwise  uniform  dense  plasmas  ? 
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Plasma  Photonic  Crystals 


Through  self-organization  of  coherent 
periodic  microstructures  of  high  plasma 
density  -  several  high  bandwidth 
reconfigurable  THz  wave  devices  can 
be  constructed: 


Directional  wave  guiding 
through  mid-band  defect 
wave  localization 


\7 


plasma  rod 
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(a) 


source 
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Directional  wave  radiation 
through  plasmon  resonances 
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continuous  wave  versus  pulsed  laser-induced  fluorescence 


Pulsed  lasers  are  transform-limited,  cannot  resolve 
distributions  in  moderately  warm  plasmas  of  relatively 
heavy  species  fe.a..  argon,  krypton,  xenon) 


Must  use  ultra-narrow  continuous 
wave  laser  sources  to  resolve 
Doppler-broadened  and  shifted 
spectral  lines  in  heavy  ions! 


*  Mazouffre  /  CNRS  used  continuous 
wave  laser  source,  was  unsuccessful 
synchronize  drifting  phases  in  the 
coherent  plasma  fluctuations 


Cappelli  /  Stanford  successfully  synchronized  continuous  wave  laser 
induced  fluorescence  to  coherent  structures  in  heavy  ion  plasmas 

•  Laser  is  chopped  at  frequency  <  coherent  fluctuation  frequency  and  individual  sample-held 
signals  are  passed  through  digital  lock-in  amplifier  to  pull  out  time-synchronized  LIF 
lineshapes  as  laser  is  scanned  in  wavelength 
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4  A  Electrospravs.  dual-mode  propulsion  (Lozano-MIT-YIP) 

Analytical  Model  reveals  distal  Contact  as  a  solution  for 
Electrochemical  Degradation  of  Emitter  Tips 


•When  a  single  polarity  is  extracted  in  the  pure  ionic  regime,  counter  ions  accumulate  in  a  double  layer  of  charge  and 
could  produce  corrosion  of  the  emitter  if  its  potential  increases  beyond  the  electrochemical  window  limit 

•Voltage  alternation  incapable  of  removing  electrochemical  degradation  at  the  tip  apex,  where  the  double  layer  potential 
grows  faster  than  its  upstream  diffusion. 


Direct  contact 

Ionic  liquid 


Wetted 

single 

emitter 


Distal  contact 

Ionic  liquid 


Ion  beam 


Extractor 


Wetted 

single 

emitter 


Distal  electrode 


HV  power 

HV  Dower 

supply 

supply 

\ 


after  80  hours  of  DC  operation 


300  pm 

'Strong  damage 
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300  pm 

intact 


Ion  beam 


Extractor 


In  addition  to  removing 
electrochemistry  from  the 
tip,  this  configuration 
allows  for  emitter 
manufacturing  using 
dielectric  materials. 


2nd  Generation: 
Thermo-Chemical  Behavior  of 

Nickel/Aluminum  Core-Shell  Particles 


•  Ignition  of  aluminum  particles  may  be  assisted  by  design  of  coating.  Examples  are 
Ti/B  and  Ni/AI  particles. 

•  What  is  the  influence  of  the  core-shell  design  on  ignition? 

Molecular  dynamics  calculations  are  performed  to  investigate  particle  designs. 


Ni  coated  Al  Particle 
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Aluminum  Core  MP 


Ni  shell  exerts 
cage-like  effect 
and  raises  MP 
of  core;  ignition 
only  after  core 
melting 


Al  coated  Ni  Particle 
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Surface  pre¬ 
melting  of  Al 
shell  implies 
lower  ignition 
temperature  of 
Al-coated  Ni 
particles. 


2.0 

Shell  Thickness,  nm 
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A  Paradigm  Shift  in  High  Pressure  Combustion  Dynamics 
Prediction  and  Control:  Analytics  and  Dynamic  Data  Driven 

Modeling  and  Validation 


•Great  Challenge  Example:  In  a  High  Pressure/Temperature,  Two-Phase,  Turbulent,  Acoustically  - 
Excited  Environment,  investigate  Amplification 

Processes: 

■Jet  Break-up 
■Atomization 
■Vaporization 
■Supercritical  States 
■Turbulence 
■Compressibility 
■Combustion 

•High  amplitude  and  high  frequency  "Acoustic  Field 

acoustic  instabilities  can  lead  to  local 

burnout  of  the  combustion  chamber  "Boundary  Interactions 

walls  and  injector  plates 

SCIENTIFIC  CHALLENGE: 

•  Modeling  and  Simulations  of  highly  complex,  nonlinear,  multi-physics,  multi-scale 
stochastic  phenomena 

•  Current  state  of  the  art  methodologies  are  not  adequate  to  address  the  challenges  in  this 
domain  new  mathematics  and  computational  methods  are  needed 
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State  of  the  Art  Modeling 


History  of  modeling  efforts  (ALL  DETERMINISTIC) 

O 

Galerkin’s  Method  /  perturbation  expansions  of  Navier  Stokes 

$ 

Unit  problems  (injector  flow,  atomization,  drop  vaporization,  spray  combustion  etc. ) 

Large  Eddy  Simulations  with  subgrid  Models 

Why  these  modeling  efforts  are  not  sufficient? 

•  Stochastic  modeling  is  essential  to  capture  the  complex  physical  phenomena  -> 
very  challenging  problem  since  current  stochastic  models  only  look  at  simplified 
problems 

•  The  problem  cannot  be  addressed  using  only  simulations,  or  experiments 


•  Use  of  experimental  data  is  essential  to  validate  models  and  codes,  but  it  should  not 
be  done  aposteriori  as  it  is  done  today 

•  Consistency  of  models  with  data  is  a  major  issue 


•  Hard  to  get  data,  either  too  much  data,  or  too  little  need  a  new  mathematical 
framework  to  bring  together  data,  experiments  and  simulations  in  a  dynamic, 
feedback  manner 
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Any  model,  by  definition,  involves  stochastic  behaviors 


(PERTURBATION  EXPANSION  EXAMPLE) 
classical  wave  equation  for  combustion  instability: 

model  uncertainties  ( the  model  has  many  approximations  and  assumptions )  including 
intrinsic  stochastic  behaviors  in  the  real  physics 

L(P')  =  dP'2  /  dt 2  -  C2dP'2/dx2  =  f(M,Q,NL ) 

X  Combustion  noise 

Parameter  uncertainty 

random  fluctuation  of  speed  of  sound  Source  term  uncertainty 

due  to  temperature  uniformity  and  turbulence-induced  random  combustion  response, 

fluctuation  acoustic  damping,  shear-layer  instability 


Parametric  Variability-  Uncertainties  at  the  Boundary  Conditions: 


PROPER  ORTHOGONAL  DECOMPOSITION  revealed  that 

Kelvin-Helmholtz  wave  motion  and  the  acoustic  waves 

must  be  accounted  as  the  boundary  conditions  for  the 

chamber  dynamics  simulation 


(^circulating  zone 


hydrodynamic  waves 
within  LOX  film 


acoustic  waves 


Kelvin-Helmholtz  instability 


Experimental  uncertainty  (dynamical  data 

sampling): 

When  the  inlet  flow  condition  fluctuates 
stochastically,  the  flame  switches  from  one 
recirculation  region  to  other  (bifurcation) 
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Compressible  LES  Equations  + 
Reduced  Basis  Model 


Stochastic  modeling  requires  many  realizations.  LES  alone  is  too  costly.  So,  a 
combined  LES-Reduced  Basis  Model  approach  fits  best. 

dp  dpuj 


Mass 

Momentum 

Energy 

Species 


+ 


dt  dxt 
dpui  d 


=  0 


H - (pUfUj  +  pSr  -  T-,  +  T-  Ss)=  0 

dt  dXj '  '  J  8  ,J  "  ’ 

“  +  J—  [u,  ( pE  +  p)  -  TvUj  +  q,  +  H’g’  +  ofs,]  =  0 

8PYk  +  8 


dt  dx; 


- pVukYk  +  Y™  +  e;y)=  wk  where  :k  =  l,Ns 


Subgrid  Stress  Flux 

TijS  =  p  ( UiUj  —  UiUj) 

Subgrid  Enthalpy  Flux 

Ht9S  =p(Eui-Eui)  +  ( UiP 

Subgrid  Viscous  Work 

=  (UjTij  -  UjTij) 

Subgrid  Species  Flux 

KT  =  P  ("'>'•■  -UiYk) 

Subgrid  Mass  Diffusion  Flux 

°'!9k  =  p  (ydkYk  -  vdkYk) 

Arrhenius  term 

k  =  Ae~EjRT 
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What  scientific  questions  does  the  proposed  BRI  address? 

Stochastic  at  Every  Scale 


•  Molecular  (atomistic)  to  micro  scales  (Angstrom-microns) 


-  Uncertainty  in  kinetics,  transport  properties 


-  Kinetics  interactions  with  fine-scale  turbulence 


•  Meso-scale  (100  micron-mm) 

-  Small-scale  turbulence-kinetics  coupling 

-  Uncertainty  in  coupling  to  large-scales 

•  Macro-scale  (cm-m) 

-  System  level  responses  and  nonlinear  feedback  effects 

-  Uncertainty  in  boundary  conditions 

•  Insufficient  data  will  lead  to  uncertainty  due  to  incomplete 
system  characterization  -  Closures  needs  to  include 
Uncertainty  Quantification  (UQ)  in  the  model 

•  Possible  UQ  Methods: 

-  Polynomial  Chaos,  Stochastic  Collocation,  Bayesian  Approaches  for  Inverse 
Modeling  and  Data  Assimilation,  Sparse  Sampling  Methods,  etc.  The  underlying 
problems  are  multiscale  and  high-dimensional.  Dealing  with  the  curse  of 
dimensionality  is  a  major  challenge. 
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Dynamic  Data  Driven  Strategic  Simulation 

&  Modeling 
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Wirz  /  UCLA-  YIP-  Center  of  Excellence  (Univ.  of  Michigan) 

Self-consistent  analytical  theory  used  to 
optimize  micro  plasma  source 


*  Critical  research  for  in-space  propulsion,  plasma-enhanced  combustion,  plasma  aerodynamics,  small  scale  sensors, 
directed  energy  devices,  plasma  processing  ... 


Goal:  Develop  efficient  and 
stable  cusp-confined  micro 
discharge 

•  New  “stream  function” 
magnetic  field  analysis 
developed  to  predict  plasma 
behavior  throughout  discharge 


Optimization  using  analytical  theory  yields  2X  increase 
in  plasma  production  over  baseline 

•  Model  trends  experimentally  validated 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited 


